3.2.12 \(\int \frac {x^3 (a+b \csc ^{-1}(c x))}{(d+e x^2)^3} \, dx\) [112]

Optimal. Leaf size=157 \[ -\frac {b c x \sqrt {-1+c^2 x^2}}{8 e \left (c^2 d+e\right ) \sqrt {c^2 x^2} \left (d+e x^2\right )}+\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {b c \left (c^2 d+2 e\right ) x \text {ArcTan}\left (\frac {\sqrt {e} \sqrt {-1+c^2 x^2}}{\sqrt {c^2 d+e}}\right )}{8 d e^{3/2} \left (c^2 d+e\right )^{3/2} \sqrt {c^2 x^2}} \]

[Out]

1/4*x^4*(a+b*arccsc(c*x))/d/(e*x^2+d)^2+1/8*b*c*(c^2*d+2*e)*x*arctan(e^(1/2)*(c^2*x^2-1)^(1/2)/(c^2*d+e)^(1/2)
)/d/e^(3/2)/(c^2*d+e)^(3/2)/(c^2*x^2)^(1/2)-1/8*b*c*x*(c^2*x^2-1)^(1/2)/e/(c^2*d+e)/(e*x^2+d)/(c^2*x^2)^(1/2)

________________________________________________________________________________________

Rubi [A]
time = 0.13, antiderivative size = 157, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 7, integrand size = 21, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.333, Rules used = {270, 5347, 12, 457, 79, 65, 211} \begin {gather*} \frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {b c x \left (c^2 d+2 e\right ) \text {ArcTan}\left (\frac {\sqrt {e} \sqrt {c^2 x^2-1}}{\sqrt {c^2 d+e}}\right )}{8 d e^{3/2} \sqrt {c^2 x^2} \left (c^2 d+e\right )^{3/2}}-\frac {b c x \sqrt {c^2 x^2-1}}{8 e \sqrt {c^2 x^2} \left (c^2 d+e\right ) \left (d+e x^2\right )} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[(x^3*(a + b*ArcCsc[c*x]))/(d + e*x^2)^3,x]

[Out]

-1/8*(b*c*x*Sqrt[-1 + c^2*x^2])/(e*(c^2*d + e)*Sqrt[c^2*x^2]*(d + e*x^2)) + (x^4*(a + b*ArcCsc[c*x]))/(4*d*(d
+ e*x^2)^2) + (b*c*(c^2*d + 2*e)*x*ArcTan[(Sqrt[e]*Sqrt[-1 + c^2*x^2])/Sqrt[c^2*d + e]])/(8*d*e^(3/2)*(c^2*d +
e)^(3/2)*Sqrt[c^2*x^2])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 65

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 79

Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.), x_Symbol] :> Simp[(-(b*e - a*f
))*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(f*(p + 1)*(c*f - d*e))), x] - Dist[(a*d*f*(n + p + 2) - b*(d*e*(n + 1
) + c*f*(p + 1)))/(f*(p + 1)*(c*f - d*e)), Int[(c + d*x)^n*(e + f*x)^(p + 1), x], x] /; FreeQ[{a, b, c, d, e,
f, n}, x] && LtQ[p, -1] && ( !LtQ[n, -1] || IntegerQ[p] ||  !(IntegerQ[n] ||  !(EqQ[e, 0] ||  !(EqQ[c, 0] || L
tQ[p, n]))))

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 270

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(c*x)^(m + 1)*((a + b*x^n)^(p + 1)/(a*
c*(m + 1))), x] /; FreeQ[{a, b, c, m, n, p}, x] && EqQ[(m + 1)/n + p + 1, 0] && NeQ[m, -1]

Rule 457

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> Dist[1/n, Subst[Int
[x^(Simplify[(m + 1)/n] - 1)*(a + b*x)^p*(c + d*x)^q, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p, q}, x] &&
NeQ[b*c - a*d, 0] && IntegerQ[Simplify[(m + 1)/n]]

Rule 5347

Int[((a_.) + ArcCsc[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_.)*((d_.) + (e_.)*(x_)^2)^(p_.), x_Symbol] :> With[{u =
IntHide[(f*x)^m*(d + e*x^2)^p, x]}, Dist[a + b*ArcCsc[c*x], u, x] + Dist[b*c*(x/Sqrt[c^2*x^2]), Int[SimplifyI
ntegrand[u/(x*Sqrt[c^2*x^2 - 1]), x], x], x]] /; FreeQ[{a, b, c, d, e, f, m, p}, x] && ((IGtQ[p, 0] &&  !(ILtQ
[(m - 1)/2, 0] && GtQ[m + 2*p + 3, 0])) || (IGtQ[(m + 1)/2, 0] &&  !(ILtQ[p, 0] && GtQ[m + 2*p + 3, 0])) || (I
LtQ[(m + 2*p + 1)/2, 0] &&  !ILtQ[(m - 1)/2, 0]))

Rubi steps

\begin {align*} \int \frac {x^3 \left (a+b \csc ^{-1}(c x)\right )}{\left (d+e x^2\right )^3} \, dx &=\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {(b c x) \int \frac {x^3}{4 d \sqrt {-1+c^2 x^2} \left (d+e x^2\right )^2} \, dx}{\sqrt {c^2 x^2}}\\ &=\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {(b c x) \int \frac {x^3}{\sqrt {-1+c^2 x^2} \left (d+e x^2\right )^2} \, dx}{4 d \sqrt {c^2 x^2}}\\ &=\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {(b c x) \text {Subst}\left (\int \frac {x}{\sqrt {-1+c^2 x} (d+e x)^2} \, dx,x,x^2\right )}{8 d \sqrt {c^2 x^2}}\\ &=-\frac {b c x \sqrt {-1+c^2 x^2}}{8 e \left (c^2 d+e\right ) \sqrt {c^2 x^2} \left (d+e x^2\right )}+\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {\left (b c \left (c^2 d+2 e\right ) x\right ) \text {Subst}\left (\int \frac {1}{\sqrt {-1+c^2 x} (d+e x)} \, dx,x,x^2\right )}{16 d e \left (c^2 d+e\right ) \sqrt {c^2 x^2}}\\ &=-\frac {b c x \sqrt {-1+c^2 x^2}}{8 e \left (c^2 d+e\right ) \sqrt {c^2 x^2} \left (d+e x^2\right )}+\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {\left (b \left (c^2 d+2 e\right ) x\right ) \text {Subst}\left (\int \frac {1}{d+\frac {e}{c^2}+\frac {e x^2}{c^2}} \, dx,x,\sqrt {-1+c^2 x^2}\right )}{8 c d e \left (c^2 d+e\right ) \sqrt {c^2 x^2}}\\ &=-\frac {b c x \sqrt {-1+c^2 x^2}}{8 e \left (c^2 d+e\right ) \sqrt {c^2 x^2} \left (d+e x^2\right )}+\frac {x^4 \left (a+b \csc ^{-1}(c x)\right )}{4 d \left (d+e x^2\right )^2}+\frac {b c \left (c^2 d+2 e\right ) x \tan ^{-1}\left (\frac {\sqrt {e} \sqrt {-1+c^2 x^2}}{\sqrt {c^2 d+e}}\right )}{8 d e^{3/2} \left (c^2 d+e\right )^{3/2} \sqrt {c^2 x^2}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [C] Result contains complex when optimal does not.
time = 0.74, size = 390, normalized size = 2.48 \begin {gather*} \frac {\frac {4 a d}{\left (d+e x^2\right )^2}-\frac {8 a}{d+e x^2}-\frac {2 b c e \sqrt {1-\frac {1}{c^2 x^2}} x}{\left (c^2 d+e\right ) \left (d+e x^2\right )}-\frac {4 b \left (d+2 e x^2\right ) \csc ^{-1}(c x)}{\left (d+e x^2\right )^2}+\frac {4 b \text {ArcSin}\left (\frac {1}{c x}\right )}{d}+\frac {b \sqrt {e} \left (c^2 d+2 e\right ) \log \left (\frac {16 d \sqrt {-c^2 d-e} e^{3/2} \left (i \sqrt {e}+c \left (c \sqrt {d}-i \sqrt {-c^2 d-e} \sqrt {1-\frac {1}{c^2 x^2}}\right ) x\right )}{b \left (c^2 d+2 e\right ) \left (\sqrt {d}+i \sqrt {e} x\right )}\right )}{d \left (-c^2 d-e\right )^{3/2}}+\frac {b \sqrt {e} \left (c^2 d+2 e\right ) \log \left (-\frac {16 d \sqrt {-c^2 d-e} e^{3/2} \left (-\sqrt {e}+c \left (-i c \sqrt {d}+\sqrt {-c^2 d-e} \sqrt {1-\frac {1}{c^2 x^2}}\right ) x\right )}{b \left (c^2 d+2 e\right ) \left (i \sqrt {d}+\sqrt {e} x\right )}\right )}{d \left (-c^2 d-e\right )^{3/2}}}{16 e^2} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[(x^3*(a + b*ArcCsc[c*x]))/(d + e*x^2)^3,x]

[Out]

((4*a*d)/(d + e*x^2)^2 - (8*a)/(d + e*x^2) - (2*b*c*e*Sqrt[1 - 1/(c^2*x^2)]*x)/((c^2*d + e)*(d + e*x^2)) - (4*
b*(d + 2*e*x^2)*ArcCsc[c*x])/(d + e*x^2)^2 + (4*b*ArcSin[1/(c*x)])/d + (b*Sqrt[e]*(c^2*d + 2*e)*Log[(16*d*Sqrt
[-(c^2*d) - e]*e^(3/2)*(I*Sqrt[e] + c*(c*Sqrt[d] - I*Sqrt[-(c^2*d) - e]*Sqrt[1 - 1/(c^2*x^2)])*x))/(b*(c^2*d +
2*e)*(Sqrt[d] + I*Sqrt[e]*x))])/(d*(-(c^2*d) - e)^(3/2)) + (b*Sqrt[e]*(c^2*d + 2*e)*Log[(-16*d*Sqrt[-(c^2*d)
- e]*e^(3/2)*(-Sqrt[e] + c*((-I)*c*Sqrt[d] + Sqrt[-(c^2*d) - e]*Sqrt[1 - 1/(c^2*x^2)])*x))/(b*(c^2*d + 2*e)*(I
*Sqrt[d] + Sqrt[e]*x))])/(d*(-(c^2*d) - e)^(3/2)))/(16*e^2)

________________________________________________________________________________________

Maple [B] Leaf count of result is larger than twice the leaf count of optimal. \(1830\) vs. \(2(135)=270\).
time = 1.48, size = 1831, normalized size = 11.66

method	result	size
derivativedivides	\(\text {Expression too large to display}\)	\(1831\)
default	\(\text {Expression too large to display}\)	\(1831\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^3*(a+b*arccsc(c*x))/(e*x^2+d)^3,x,method=_RETURNVERBOSE)

[Out]

1/c^4*(a*c^6*(-1/2/e^2/(c^2*e*x^2+c^2*d)+1/4*d*c^2/e^2/(c^2*e*x^2+c^2*d)^2)-1/2*b*c^6*arccsc(c*x)/e^2/(c^2*e*x
^2+c^2*d)+1/4*b*c^8*arccsc(c*x)*d/e^2/(c^2*e*x^2+c^2*d)^2-1/4*b*c^7*(c^2*x^2-1)^(1/2)/e/((c^2*x^2-1)/c^2/x^2)^
(1/2)/x*d/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*arctan(1/(c^2*x^2-1)^(1/2))-1/4*b*c^7*(
c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*ar
ctan(1/(c^2*x^2-1)^(1/2))+1/16*b*c^7*(c^2*x^2-1)^(1/2)/e/((c^2*x^2-1)/c^2/x^2)^(1/2)/x*d/(-(c^2*d+e)/e)^(1/2)/
(c^2*d+e)/(-e*c*x+(-c^2*d*e)^(1/2))/(e*c*x+(-c^2*d*e)^(1/2))*ln(2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e-(-
c^2*d*e)^(1/2)*c*x-e)/(e*c*x+(-c^2*d*e)^(1/2)))+1/16*b*c^7*(c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/(-(
c^2*d+e)/e)^(1/2)/(c^2*d+e)/(-e*c*x+(-c^2*d*e)^(1/2))/(e*c*x+(-c^2*d*e)^(1/2))*ln(2*((-(c^2*d+e)/e)^(1/2)*(c^2
*x^2-1)^(1/2)*e-(-c^2*d*e)^(1/2)*c*x-e)/(e*c*x+(-c^2*d*e)^(1/2)))+1/16*b*c^7*(c^2*x^2-1)^(1/2)/e/((c^2*x^2-1)/
c^2/x^2)^(1/2)/x*d/(-(c^2*d+e)/e)^(1/2)/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*ln(-2*((-
(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e+(-c^2*d*e)^(1/2)*c*x-e)/(-e*c*x+(-c^2*d*e)^(1/2)))+1/16*b*c^7*(c^2*x^2-
1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/(-(c^2*d+e)/e)^(1/2)/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d
*e)^(1/2))*ln(-2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e+(-c^2*d*e)^(1/2)*c*x-e)/(-e*c*x+(-c^2*d*e)^(1/2)))+
1/8*b*c^5*(c^2*x^2-1)/((c^2*x^2-1)/c^2/x^2)^(1/2)/x/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2
))-1/4*b*c^5*(c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)/x/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*
d*e)^(1/2))*arctan(1/(c^2*x^2-1)^(1/2))-1/4*b*c^5*(c^2*x^2-1)^(1/2)*e/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/d/(c^2*d+e
)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*arctan(1/(c^2*x^2-1)^(1/2))+1/8*b*c^5*(c^2*x^2-1)^(1/2)/(
(c^2*x^2-1)/c^2/x^2)^(1/2)/x/(-(c^2*d+e)/e)^(1/2)/(c^2*d+e)/(-e*c*x+(-c^2*d*e)^(1/2))/(e*c*x+(-c^2*d*e)^(1/2))
*ln(2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e-(-c^2*d*e)^(1/2)*c*x-e)/(e*c*x+(-c^2*d*e)^(1/2)))+1/8*b*c^5*(c
^2*x^2-1)^(1/2)*e/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/d/(-(c^2*d+e)/e)^(1/2)/(c^2*d+e)/(-e*c*x+(-c^2*d*e)^(1/2))/(e*
c*x+(-c^2*d*e)^(1/2))*ln(2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e-(-c^2*d*e)^(1/2)*c*x-e)/(e*c*x+(-c^2*d*e)
^(1/2)))+1/8*b*c^5*(c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)/x/(-(c^2*d+e)/e)^(1/2)/(c^2*d+e)/(e*c*x+(-c^2
*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*ln(-2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e+(-c^2*d*e)^(1/2)*c*x-e)
/(-e*c*x+(-c^2*d*e)^(1/2)))+1/8*b*c^5*(c^2*x^2-1)^(1/2)*e/((c^2*x^2-1)/c^2/x^2)^(1/2)*x/d/(-(c^2*d+e)/e)^(1/2)
/(c^2*d+e)/(e*c*x+(-c^2*d*e)^(1/2))/(-e*c*x+(-c^2*d*e)^(1/2))*ln(-2*((-(c^2*d+e)/e)^(1/2)*(c^2*x^2-1)^(1/2)*e+
(-c^2*d*e)^(1/2)*c*x-e)/(-e*c*x+(-c^2*d*e)^(1/2))))

________________________________________________________________________________________

Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(a+b*arccsc(c*x))/(e*x^2+d)^3,x, algorithm="maxima")

[Out]

-1/4*(2*x^2*e + d)*a/(x^4*e^4 + 2*d*x^2*e^3 + d^2*e^2) - 1/4*(2*x^2*arctan2(1, sqrt(c*x + 1)*sqrt(c*x - 1))*e
+ d*arctan2(1, sqrt(c*x + 1)*sqrt(c*x - 1)) + 4*(x^4*e^4 + 2*d*x^2*e^3 + d^2*e^2)*integrate(1/4*(2*c^2*x^3*e +
c^2*d*x)*e^(1/2*log(c*x + 1) + 1/2*log(c*x - 1))/(c^2*x^6*e^4 + (2*c^2*d*e^3 - e^4)*x^4 + (c^2*d^2*e^2 - 2*d*
e^3)*x^2 - d^2*e^2 + (c^2*x^6*e^4 + (2*c^2*d*e^3 - e^4)*x^4 + (c^2*d^2*e^2 - 2*d*e^3)*x^2 - d^2*e^2)*e^(log(c*
x + 1) + log(c*x - 1))), x))*b/(x^4*e^4 + 2*d*x^2*e^3 + d^2*e^2)

________________________________________________________________________________________

Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 496 vs. \(2 (138) = 276\).
time = 0.60, size = 1021, normalized size = 6.50 \begin {gather*} \left [-\frac {4 \, a c^{4} d^{4} + 8 \, a d x^{2} e^{3} + {\left (b c^{2} d^{3} + 2 \, b x^{4} e^{3} + {\left (b c^{2} d x^{4} + 4 \, b d x^{2}\right )} e^{2} + 2 \, {\left (b c^{2} d^{2} x^{2} + b d^{2}\right )} e\right )} \sqrt {-c^{2} d e - e^{2}} \log \left (-\frac {c^{2} d - {\left (c^{2} x^{2} - 2\right )} e + 2 \, \sqrt {c^{2} x^{2} - 1} \sqrt {-c^{2} d e - e^{2}}}{x^{2} e + d}\right ) + 4 \, {\left (b c^{4} d^{4} + 2 \, b d x^{2} e^{3} + {\left (4 \, b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2} + 2 \, {\left (b c^{4} d^{3} x^{2} + b c^{2} d^{3}\right )} e\right )} \operatorname {arccsc}\left (c x\right ) + 8 \, {\left (b c^{4} d^{4} + b x^{4} e^{4} + 2 \, {\left (b c^{2} d x^{4} + b d x^{2}\right )} e^{3} + {\left (b c^{4} d^{2} x^{4} + 4 \, b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2} + 2 \, {\left (b c^{4} d^{3} x^{2} + b c^{2} d^{3}\right )} e\right )} \arctan \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) + 4 \, {\left (4 \, a c^{2} d^{2} x^{2} + a d^{2}\right )} e^{2} + 8 \, {\left (a c^{4} d^{3} x^{2} + a c^{2} d^{3}\right )} e + 2 \, {\left (b c^{2} d^{3} e + b d x^{2} e^{3} + {\left (b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2}\right )} \sqrt {c^{2} x^{2} - 1}}{16 \, {\left (c^{4} d^{5} e^{2} + d x^{4} e^{6} + 2 \, {\left (c^{2} d^{2} x^{4} + d^{2} x^{2}\right )} e^{5} + {\left (c^{4} d^{3} x^{4} + 4 \, c^{2} d^{3} x^{2} + d^{3}\right )} e^{4} + 2 \, {\left (c^{4} d^{4} x^{2} + c^{2} d^{4}\right )} e^{3}\right )}}, -\frac {2 \, a c^{4} d^{4} + 4 \, a d x^{2} e^{3} - {\left (b c^{2} d^{3} + 2 \, b x^{4} e^{3} + {\left (b c^{2} d x^{4} + 4 \, b d x^{2}\right )} e^{2} + 2 \, {\left (b c^{2} d^{2} x^{2} + b d^{2}\right )} e\right )} \sqrt {c^{2} d e + e^{2}} \arctan \left (\frac {\sqrt {c^{2} x^{2} - 1} \sqrt {c^{2} d e + e^{2}}}{c^{2} d + e}\right ) + 2 \, {\left (b c^{4} d^{4} + 2 \, b d x^{2} e^{3} + {\left (4 \, b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2} + 2 \, {\left (b c^{4} d^{3} x^{2} + b c^{2} d^{3}\right )} e\right )} \operatorname {arccsc}\left (c x\right ) + 4 \, {\left (b c^{4} d^{4} + b x^{4} e^{4} + 2 \, {\left (b c^{2} d x^{4} + b d x^{2}\right )} e^{3} + {\left (b c^{4} d^{2} x^{4} + 4 \, b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2} + 2 \, {\left (b c^{4} d^{3} x^{2} + b c^{2} d^{3}\right )} e\right )} \arctan \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) + 2 \, {\left (4 \, a c^{2} d^{2} x^{2} + a d^{2}\right )} e^{2} + 4 \, {\left (a c^{4} d^{3} x^{2} + a c^{2} d^{3}\right )} e + {\left (b c^{2} d^{3} e + b d x^{2} e^{3} + {\left (b c^{2} d^{2} x^{2} + b d^{2}\right )} e^{2}\right )} \sqrt {c^{2} x^{2} - 1}}{8 \, {\left (c^{4} d^{5} e^{2} + d x^{4} e^{6} + 2 \, {\left (c^{2} d^{2} x^{4} + d^{2} x^{2}\right )} e^{5} + {\left (c^{4} d^{3} x^{4} + 4 \, c^{2} d^{3} x^{2} + d^{3}\right )} e^{4} + 2 \, {\left (c^{4} d^{4} x^{2} + c^{2} d^{4}\right )} e^{3}\right )}}\right ] \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(a+b*arccsc(c*x))/(e*x^2+d)^3,x, algorithm="fricas")

[Out]

[-1/16*(4*a*c^4*d^4 + 8*a*d*x^2*e^3 + (b*c^2*d^3 + 2*b*x^4*e^3 + (b*c^2*d*x^4 + 4*b*d*x^2)*e^2 + 2*(b*c^2*d^2*
x^2 + b*d^2)*e)*sqrt(-c^2*d*e - e^2)*log(-(c^2*d - (c^2*x^2 - 2)*e + 2*sqrt(c^2*x^2 - 1)*sqrt(-c^2*d*e - e^2))
/(x^2*e + d)) + 4*(b*c^4*d^4 + 2*b*d*x^2*e^3 + (4*b*c^2*d^2*x^2 + b*d^2)*e^2 + 2*(b*c^4*d^3*x^2 + b*c^2*d^3)*e
)*arccsc(c*x) + 8*(b*c^4*d^4 + b*x^4*e^4 + 2*(b*c^2*d*x^4 + b*d*x^2)*e^3 + (b*c^4*d^2*x^4 + 4*b*c^2*d^2*x^2 +
b*d^2)*e^2 + 2*(b*c^4*d^3*x^2 + b*c^2*d^3)*e)*arctan(-c*x + sqrt(c^2*x^2 - 1)) + 4*(4*a*c^2*d^2*x^2 + a*d^2)*e
^2 + 8*(a*c^4*d^3*x^2 + a*c^2*d^3)*e + 2*(b*c^2*d^3*e + b*d*x^2*e^3 + (b*c^2*d^2*x^2 + b*d^2)*e^2)*sqrt(c^2*x^
2 - 1))/(c^4*d^5*e^2 + d*x^4*e^6 + 2*(c^2*d^2*x^4 + d^2*x^2)*e^5 + (c^4*d^3*x^4 + 4*c^2*d^3*x^2 + d^3)*e^4 + 2
*(c^4*d^4*x^2 + c^2*d^4)*e^3), -1/8*(2*a*c^4*d^4 + 4*a*d*x^2*e^3 - (b*c^2*d^3 + 2*b*x^4*e^3 + (b*c^2*d*x^4 + 4
*b*d*x^2)*e^2 + 2*(b*c^2*d^2*x^2 + b*d^2)*e)*sqrt(c^2*d*e + e^2)*arctan(sqrt(c^2*x^2 - 1)*sqrt(c^2*d*e + e^2)/
(c^2*d + e)) + 2*(b*c^4*d^4 + 2*b*d*x^2*e^3 + (4*b*c^2*d^2*x^2 + b*d^2)*e^2 + 2*(b*c^4*d^3*x^2 + b*c^2*d^3)*e)
*arccsc(c*x) + 4*(b*c^4*d^4 + b*x^4*e^4 + 2*(b*c^2*d*x^4 + b*d*x^2)*e^3 + (b*c^4*d^2*x^4 + 4*b*c^2*d^2*x^2 + b
*d^2)*e^2 + 2*(b*c^4*d^3*x^2 + b*c^2*d^3)*e)*arctan(-c*x + sqrt(c^2*x^2 - 1)) + 2*(4*a*c^2*d^2*x^2 + a*d^2)*e^
2 + 4*(a*c^4*d^3*x^2 + a*c^2*d^3)*e + (b*c^2*d^3*e + b*d*x^2*e^3 + (b*c^2*d^2*x^2 + b*d^2)*e^2)*sqrt(c^2*x^2 -
1))/(c^4*d^5*e^2 + d*x^4*e^6 + 2*(c^2*d^2*x^4 + d^2*x^2)*e^5 + (c^4*d^3*x^4 + 4*c^2*d^3*x^2 + d^3)*e^4 + 2*(c
^4*d^4*x^2 + c^2*d^4)*e^3)]

________________________________________________________________________________________

Sympy [F(-1)] Timed out
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**3*(a+b*acsc(c*x))/(e*x**2+d)**3,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F(-2)]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: RuntimeError} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3*(a+b*arccsc(c*x))/(e*x^2+d)^3,x, algorithm="giac")

[Out]

Exception raised: RuntimeError >> An error occurred running a Giac command:INPUT:sage2OUTPUT:Warning, integrat
ion of abs or sign assumes constant sign by intervals (correct if the argument is real):Check [abs(sageVARx)]s
ym2poly/r2sym(

________________________________________________________________________________________

Mupad [F]
time = 0.00, size = -1, normalized size = -0.01 \begin {gather*} \int \frac {x^3\,\left (a+b\,\mathrm {asin}\left (\frac {1}{c\,x}\right )\right )}{{\left (e\,x^2+d\right )}^3} \,d x \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((x^3*(a + b*asin(1/(c*x))))/(d + e*x^2)^3,x)

[Out]

int((x^3*(a + b*asin(1/(c*x))))/(d + e*x^2)^3, x)

________________________________________________________________________________________